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ABSTRACT 



Aims. To study the mass distribution of galaxy clusters up to their Virial radius, CL0016+16 seems to be a good 
candidate, since it is a bright massive cluster, previously considered as being dynamically relaxed. 
Methods. Using XMM-Newton observations of CL0016+16, we performed a careful X-ray background analysis, and, we 
detected convincingly its X-ray emission up to R200 ■ We then studied its dynamical state with a detailed 2D temperature 
and surface brightness analysis of the inner part of the cluster. 

Results. Using the assumption of both spherical symmetry and hydrostatic equilibrium (HE) we can determine the 
main cluster parameters; total mass, temperature profile, surface brightness profile and /3-parameter. We also build a 
temperature map which clearly exhibits departure from spherical symmetry in the centre. To estimate the influence of 
these perturbations onto our total mass estimate, we also compute the total mass in the framework of the HE approach, 
but this time with various temperature profiles obtained in different directions. These various total mass estimates are 
consistent with each other. The temperature perturbations are clear signatures of ongoing merger activity. We also find 
significant residuals after subtracting the emissivity map by a 2D /3-model fit. We conclude that, although CL0016-I-16 
shows clear signs of merger activity and departure from spherical symmetry in the centre, its X-ray emissivity can be 
detected up to -R200 and the corresponding mass M200 can be computed directly. It is therefore a good candidate to 
study cosmological scaling laws as predicted by the theory. 

Key words. galaxies:cluster:individual: CL0016-I-16: observation-X-rays 



1. Introduction 

In the hierarchical scenario of structure formation, clusters 
of galaxies are the largest and youngest virialized objects in 
the Universe. This makes them ideal targets for cosmologi- 
cal studies. Clusters of galaxies are self-similar in shape, and 
the cluster population obeys scaling laws for various physi- 
cal properties : total mass, temperature and luminosity. The 
evolution of these properties with redshift gives also com- 
plementary information that shed light on cluster physics. 
X-ray observations of galaxy clusters allow us to study the 
hot intra cluster medium (ICM) , which is the main baryon 
reservoir in galaxy clusters. In the era of XMM-Newton and 
Chandra observations, we can obtain detailed information 
on the density and temperature distribution of the ICM, 
and study their internal structure with unprecedented ac- 
curacy. 

This paper reports a detailed study of CL0016-fl6, 
a very massive, luminous and distant (z = 0.54) 
galaxy cluster. CL0016-I-16 is one of the most ex- 
tensively studied cluster of galaxies in d ifferent wave- 
lengt hs, in p articular in the op ti cal ()Tanaka et aTl . 
|2005[) . radio (jOiovannini fc Ferg tti'. '200 QD ban d s an d 
for Sunyaev-Zel'd o vich effect (|Birkinshaw et all . 119811) . 



The question about the dynamical state of this clus- 
ter is still open. Previous X-r ay analysis conclud e d tha t 
the cluster might be relaxed (jKotov fc Vikhlininl . l2005l ). 
But we know that t his c luster has a radio halo from 
iGiovannini fc Ferettil (|2000D . which can be a signature of 
merger activity. Since the total mass is computed using the 
hydrostatic equilibrium approach, we need to have a better 
idea on the dynamical state of the cluster. 

In this paper, we present a detailed spectro-imaging 
study of the galaxy cluster CL0016-I-16. We first assumed 
that the cluster is relaxed and spherically symmetric. We 
obtained high quality surface brightness profiles using dif- 
ferent backgrounds subtraction methods and demonstrated 
that we can detect the cluster emission up to i?200- This 
specific radius was first used by theoreticians to define dark 
matter halos, for which the mean density is 200 times the 
critical density of the universe. N body and hydrodynamics 
numerical simulation have shown that within this radius, 
the gas and dark matter halo can be reasonably considere d 
as being in dynamical equilibrium (|Cole fc Lacevl . |1996[ ). 
To compute R200, we use the standard definition: 



An 



3 

200 



(|Bonamente et al.l 120061). It h a s bee n also studied us- 



ing wea k lensing [ Cl owe et a"D ('2000!) and X-rays from 
ROSAT (lNeumann~fc Bohringcr, 1997 ) and XMM-Newto n 
(jKotov fc Vikhlininl . 120051) . (jWorraU fc Birkinshawl [200I . 
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where pcrit = SHq/SttG is the critical density, our cos- 
mological parameters — 0.3, flA — 0.7 and Hq — 70 
km s~^Mpc~"'^. Most theoretical predictions (mass function, 
density profiles, scaling laws) are computed within this ra- 
dius, while most X-ray observations are limited to much 
inner regions, such as i?iooo- The purpose of this paper is 
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to compute the observed properties of a large X-ray clus- 
ter within i?200 ■ Using the hydrostatic equilibrium and first 
guesses of the /3-model parameters and mean tem perature 
of 8.9 keV taken from iKotov fc VikhlininI (|2005[) . we esti- 
mated using Eqdithat i?200 = 1-92 Mpc or 5' for z=0.54. 
We, then verify that our measurement of temperature at 
i?200 is very similar to the used value. In a second step, we 
performed a detailed spectral analysis of CL0016-I-16, and 
deduce its average temperature, its temperature profile and 
a high quality temperature map. To estimate the influence 
of temperature variations on the total mass estimate, we 
calculated the mass profile using the hydrostatic equilib- 
rium equation and different temperature profiles in different 
directions. Finally, using the ACDM cosmological model, 
we have tested how CL0016-I-16 fits into cosmological self- 
similarity theory. All errors on the cluster parameters were 
obtained at the 68% confidence level. 



2. Data analysis 

We studied the physical parameters and ICM dynamics of 
CL0016-I-16 up to i?200- The cluster emission is weak in the 
external region, so the astrophysical background can play 
the main role in the outer region, that is why the treatment 
of background is very important. To obtain the best result 
we used three different backgrounds for subtraction in the 
data analysis. These were the background of A. Read with- 
out "flare" rejection, the background of J. Nevalainen with 
"flare" rejection and the modelling of background using the 
observation data. 

The meth o d of double background subtraction from 
lArnaud et al.l (|2002f) was used. For analysis we used the 
XMM- Newton data from EPIC cameras (MOS 1,2 and 
pn) and the XMM-Newton Science Analysis System (SAS) 
for data reduction. In the MOS 1,2 data set we took into 
account event patterns to 12 and in pn data - patterns 
to 4, flag = 0. 

The sky coordinates of background observation in the 
event files were modified using the aspect solution of 
CL0016-I-16 observation. 

From count rate of the observation data we detected 
and excluded the periods of "flare". Using the count rate 
of the observation in the high energy bands (10-12 keV) 
and exposure time of the observation we normalized the 
background. 

The effective area of the XMM-Newton mirrors is a 
function of the off-axis angle and energy of the photons. One 
could manage it thro ugh a weight func t ion d irectly com- 
puted for each event (jMaierowicz et al] . l2002f ) or through 
exposure map. The vignetting is a geometrical effect, both 
approaches lead to similar results. All along of this work 
we have used the weight method. 

We excluded all detectable point sources from the data 
observation in our spectral and spatial analysis. The sources 
were detected in the 0.3-4.5 keV energy bands. Detected 
point sources were masked with circles of 70% point spread 
function power radii. 

2.1. Double subtraction 

After cleaning the flare events, the XMM background is 
dominated by the Cosmic X-ray background (CXB) and 
non X-ray background (NXB) induced by high energy parti- 



cle s. In our analy s is we u sed the double-subtraction method 
bv lArnaud et alj ()2002 f) for background subtraction. 

The background subtraction for profile and spectrum 
consists of two steps. Firstly we subtract the normalized 
blank field obtained using the same spatial and energy se- 
lection (NXB component) and then we subtract the residual 
components using the data in the outer part of the cluster 
emission (CXB component). 

2.2. Background subtraction (blank field of A. Read) 

The background of A. Read is a blank fiel d without the 
any "fiare" rejection (jRead fc Ponmanl . l2003f) . We excluded 
the "flare" period using the similar method for the ob- 
servation and background. We suggested that the residual 
"fiare" background could play a role in each energy band. 
That is why for the observation and background data we 
performed the "fiare" re jection in the standard selection 
(jMaierowicz et all |2002[ ) and for each energy band (10-12 
keV, 0.3-12 keV, 2-5 keV, 0.3-2 keV). To determine the 
best limit of "flare" background we used the fight curve, 
we compared the histograms obtained from the images in 
the external regions. Using different energy bands for the 
"flare" rejection we obtained the same results and the same 
exposure times. 



2.3. Background subtraction (blank field of J. Nevalainen) 

To obtain the best result from the image analysis, to detect 
the cluster emission up to ii!200 with the XMM-Newton data 
we also performe d the image analysis wit h the background 
of J. Nevalainen (jNevalainen et aUl2005f) . The background 
data of J. Nevalainen is cleaned for "ffares", has better 
statistic, less sources, no artefacts in the centre of FOV. 
In this analysis for observation we excluded the "flare" re- 
jection in the standard selection using the Poisson filter the 
fiares were detected as > Str deviation from the mean. 



3. Image analysis 

3.1. Surface brightness profile 

The surface brightness profile was fi tted with the /3-model 
(jCavaliere and Fusco-Femianol Il976i) in which the surface 
brightness S{r) is defined as: 



S{r) ^So[l + ^ 



-3/3+0.5 



(2) 



where Sq is the central intensity, Vf. is the core radius, (3 
is the slope parameter. The best fit /3-model was used to 
compute the total mass. The /3-model is a model that allows 
us to project easily the surface brightness of the ICM which 
emits via thermal bremsstrahlung into the following gas 
density profile: 



n{r) = no [1 + ^ 



2 \ -3/3/2 



(3) 



3.1.1. Surface brightness profile with blank field of A. Read 

To obtain the surface brightness profile using the back- 
ground of A. Read we created the images in the 0.3-4.5 
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Fig. 1. The surface brightness profiles combined M0S1,M0S2 and pn cameras obtained using different backgrounds a) 
the background of A. Read, b) the background of J. Nevalainen, c) the background model from observation data. Solid 
line is the best fit /3-model. 



keV and calculated the corresponding exposure maps tak- 
ing into account the detector geometry for the observation 
and background data for three cameras of XMM-Newton. 
The radial surface brightness profiles were created from 
these images and three profiles were summed. Significant 
cluster emission was detected up to i?200 with the limit of 
detection 3a. The /3-fits of cluster profile are given in the 
Table [T] and the convolved best fit /3-model is plotted in 
FigM- 

3.1.2. Surface brightness profile with blank field of J. 
Nevalainen 

To obtain surface brightness profile using background of J. 
Nevalainen we extracted the surface brightness profile of 
the cluster in the 0.3-4.5 keV energy band. This band was 
chosen to optimise the signal-to-noise ratio. We binned the 
photons into concentric annuli with a size of 1.65" centred 
on the maximum of the X-ray emission for each camera. 
The three profiles were then summed. The resulting sur- 
face brightness profile S'o is shown in Fig. [TJd. The cluster 
emission is clearly detected up to i?200> with limit of the 
detection 3cr. We fitted S'o with /3-model, without the first 
three points and including the PSF deconvolution (see the 
results in Table [1]). The same results of the /3-model fit pa- 
rameters were obtained using the all points for the fits, /3 
= (0.73 ± 0.02) and = (0.27 ± 0.02) Mpc. This result is 
in a g ood agreement with the result of iKotov fc VikhlininI 
(I2005f) . 

3.1.3. Surface brightness profile using the modelling of 
background from observation data 

The astrophysical background plays a very important role 
in the data analysis of the outskirts of galaxy clusters so 
we assumed that would be better to use the observation 
data for the background modelling. Using observation data 
we modelled the CXB and NXB components of the back- 
ground. The CXB component of backgrou nd is vignetted 
by X- optics, but the NXB component is not (jArnaud et al.l . 
I2OO2I) . In this analysis we did the following: selection of the 
"flare" events in the observation data, the surface bright- 
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Fig. 2. The combined (M0S1,M0S2, pn) surface bright- 
ness profiles in the external region obtained using different 
backgrounds a) the background of A. Read, b) the back- 
ground of J. Nevalainen, c) the background modelling from 
corresponding observation data. Solid line is the best-fit /3- 
model. 



ness profile creation from observation data, using the vi- 
gnetting function we searched the CXB and NXB com- 
ponents of the background by minimising in the outer 
regions. Finally we subtracted the CXB and NXB compo- 
nents from the observation. We obtained cluster emission 
up to i?200- After rebinning the surface brightness profile 
was fitted with /3-model, see FigHJ;. The fit of the cluster 
profile is given in Table [1] 

3.2. Summary of image analysis 

In our image analysis, the same results were obtained using 
different backgrounds for subtraction. The surface bright- 
ness profile of CL0016+16 was detected up to R200, this is 
important for studying the cluster physics in the outskirts. 

In all samples the /3-model is in a good agreement with 
the surface brightness profiles using the method of double 
background subtraction. The different results of the fitting 
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Table 1. The comparison of the /3-model fitting results for the galaxy cluster CL0016+16. 





p-nt 


t^U.Ui-Z.O Wipe j 




/3-fit by region extern 


[KJ.o-Z.O MpCj 




LlcLLcL bULliCc 


P 




A.reduced 


P 




\reduced 


Kotov and Vikhlinin 2005 


0.76 ± 0.01 


0.27 ± 0.01 


1.38 








Worrall and Birkinshaw 2003 


0.70 ± 0.01 


0.23 ± 0.01 


1.35 








Neumann and Bohringer 1996 


0.68 ± 0.3 


0.28 ± 0.14 


1.02 








Using background by A. Read 


0.77 ± 0.01 


0.31 ± 0.01 


1.56 


0.80 ± 0.02 


0.34 ± 0.02 


1.35 


Using background by J.Nevalainen 


0.76 ± 0.01 


0.29 ± 0.01 


1.15 


0.72 ± 0.01 


0.29 =b 0.01 


1.02 


Using the modelling 














of background from observation 


0.76 ± 0.01 


0.38 ± 0.01 


1.57 


0.81 ± 0.02 


0.38 =b 0.02 


1.21 



/3-niodel for CL0016+16 are shown in Table[Tl a good agree- 
ment with other authors was obtained. 

We obtained worse reduced using the background 
of A. Read and the modelling of background. It may be 
because background of A. Read has artefacts in the centre 
of FOV so the surface brightness profile of the cluster does 
not fit well in the centre with /3-model. Also in this case the 
variation of reduced depends on the second background 
subtraction. 

We established that the /3-fit had problems in fitting 
the core radius, so we fitted the /3-model with the surface 
brightness profile in the external region and obtained the 
better values of reduced x^- Also, using the template of 
simulated cluster we tested the possibility of detection up 
to -R200 with XMM-Newton data (see Appendix \X\ ) . 

With XMM-Newton data we detect cluster emission of 
CL0016-I-16 up to i?200 then it was possible to determine 
the cluster total mass and physical parameters more pre- 
cisely, to study physics near the virial radius and to test 
self similarity theory. For our spectral analysis we decided 
to use the background of Nevalainen, because using this 
background we obtained a better result with the /3-model 
fit. 



4. Spectral analysis 

We needed to study the ICM dynamics of CL0016-I-16 and 
to calculate more precisely the total mass profile up to i?200 
using the equilibrium approach. 

From the spherical brightness emission, CL0016-I-16 
looks as a relaxed cluster without any cool core. On the 
other hand C L0016-I-16 has the companio n cluster RX 
J001 8.3-H618 (IWorrall k Birkinshawl I2003D and halo ra- 
dio Jg iovannini &: Ferettil . |2000| ). It is very important to 
check the equilibrium of this system. 

To test the equilibrium we performed a detailed spectral 
analysis. The mean temperature, temperature profile, tem- 
perature map, temperature in the regions and the temper- 
ature profiles in different directions were obtained. In our 
spectral analysis we used three cameras of XMM-Newton, 
the background from J. Nevalain en and the me t hod of dou- 
ble background subtraction by lArnaud et al.1 (|2002f) . We 
fitted the spectra with XSPEC using the redshifted APEC 
plasma emission model with the absorption Nh — 4 • 10^° 
cm~^ and free abundances. In our research we are more 



interested in temperature variations, at same time we fit- 
ted temperature with free and fixed abundance (0.18 solar 
unit, obtained from the best fit of the mean temperature in 
6.4' radius). For each spectrum similar temperatures were 
obtained in both case but sometimes we cannot determine 
reliably the abundances due to weak statistics. 



4.1. Mean temperature 

The overall MOSl, M0S2 and pn spectra extracted from 
the event file are shown in Fig. [3] The spectra were cor- 
rected for the vignetting and the background. The integra- 
tion region for the cluster was restricted to 5' (i?20o): this 
region was chosen to test the self similarity theory. The 
temperature values were estimated for all EPIC cameras 
MOSl, M0S2 and pn. The best fit gives kT = 8.81 ± 0.35 
keV and an abundance of 0.18 ± 0.05. The reduced of 
0.93. We also estimated the mean temperature in 6.4' ra- 
dius, the best fit gives kT = 8.83 ± 0.36 keV with reduced 
of 0.93 and an abundance of of 0.19 ± 0.05 solar units. 




channel energy (keV) 

Fig. 3. CL0016+16 extracted spectra, in the region 5', 
from the XMM-Newton data. Black, red, green points are 
the data from the corresponding EPIC/ MOSl, M0S2 and 
pn data. The solid lines is the best fit of the isothermal 
model with kT = 8.81 keV and abundance of 0.18 solar 
units. 



4.2. Temperature profile 

For a more exact determination of the total mass profile it 
is better to use the temperature profile. We supposed that 
CL0016-I-16 was a relaxed cluster and assumed that the 
temperature structure of this cluster was spherically sym- 
metric. The spectra were extracted in the five concentric 
annuli centred on the cluster X-ray emission peak and we 
fitted the data as described above. The sizes of the annuli 
were chosen to optimise the signal-to-noise in each annulus. 
Note, we fitted spectra for temperature profile with frozen 
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Table 2. The spectral fits results of the temperature profile with frozen and free abundances for all EPIC cameras, the 
abundance fits, the temperature for MOSl and M0S2 cameras, pn camera. 



annulus T{keV)/xted 


nkeV)/xied Z ( Zq) 


T[keV)/xied 


nkeV)/xied 


Z = 0.182'q 


free Z 


free Z 


free Z 


''i-''2(') niosl+nios2+pn 


mosl-|-nios2+pn mosl+mos2+pn 


mosl+mos2 


pn 



0.0-0.5 10.63lo j7/0.92 lQ.27tl%/Q-9l 0.42 ±0.11 11.45^0.95/0.88 9.58to ji/0.93 

0.5-1.0 9.3llE;:g/1.01 9.34^0.43/1.01 0.15 ±0.06 9.80t^;|?/1.02 9.05t^;g6/1.00 

1.0-2.0 9.17+0-51/0.92 9.33+054/0.91 0.11 ±0.07 10.23+j-8g/0.94 8.77+0 63 /0.89 

2.0-4.0 6.45tE;:^|/1.03 6.32t°;^5/1.03 0.38 ±0.17 8.06l:i;|^/1.03 5.19t;|! j^/1.02 

4.0-6.4 4.14j;^'^*/0.87 4.07l:i'3j/0.85 ... 1.99to 55/O.89 4.98^2' gVo.82 



and free abundances we obtained a similar values of the 
temperature and error bar (see Table [2]). 

We obtained the total temperature profile up to i?2oo 
(see Fig|4]). The temperat ure profile is in goo d agreement 
with spectral fit results of lKotov fc VikhlininI |2005.) . The 
solid line profile is that f rom our analysis and the dotted 
line profile is the profile bv lKotov fc VikhlininI (|2005( ) . Note, 
that temperature was detected up to i?200 and obtained 
points are in a go od agreement with the fun ction of tem- 
perature profile bv lKotov fc VikhlininI (|200£f) . 




I I LJ 

0.1 1.0 

R, Mpc 

Fig. 4. The solid points are the total temperatures best 
fit, the solid line is the fitti ng of the temperat u re pro - 
file with the function from iKotov fc VikhlininI (|2005l ). 
the thin points ar e the total temperatures best fit of 
iKotov fc VikhlininI (|2005l ) a nd dotted line are the obta ined 
total temperatures profile of lKotov fc VikhlininI (|2005( ). 



4.3. Temperature map 

To understand the ICM dynamics of CL0016+16 it is neces- 
sary to obtain a temperature map, which is the best accessi- 
ble measurable indicator of a system equilibrium. The tem- 
perature map was obtained by X-ray wavelet spectral map- 
ping algorithm (XWSM). The spatial temperature varia- 
tions are coded at different scales in the wavelet space using 



the Haar wavelet and denoised by thresholding the wavelet 
coe fficients. For a comp lete description of the algorithm, 
see iBourdin et all (|2004D . 

The CL0016±16 temperature map was computed using 
three cameras of XMM-Newton. This algorithm fitted the 
temperature with the APEC plasma emission model. The 
resulting temperature map shown in Fig|5] was obtained 
from a wavelet analysis performed on 5 scales corresponding 
to a structure of a minimum size of 45" using the mean 
temperature of 9 keV and with the same absorption and 
abundance. 

The temperature map was obtained up to 2.5' (1 Mpc) 
with a 68% confidence level for all detected structures. The 
overall appearance of the temperature structure is strongly 
asymmetric. The maximum of the temperature is 12.52 keV 
and in the cold region of the cluster it is 8.7 keV. We ob- 
served the maximum of temperature on the south-west from 
the cluster emission centre, and minimum of temperature 
on north-east from cluster emission centre. Fig.[5]shows the 
hot regions to south-east and to north-west and two cold 
regions in other directions. Note, that the maximum of the 
temperature is not superimposed onto the maximum emis- 
sion. The peak of the temperature is displaced by the 0.5' 
(190 kpc) from the cluster emission peak. 




Fig. 5. a) the X-ray emission from EPIC/MOSl, M0S2 
fcpn in the energy band 0.3-4.5 keV. b) the temper- 
ature map obtained by wavelet algorithm described in 
IBourdin et aLl (|2004D . 



4.4. Temperature in selected regions 

To check the temperature map obtained with XWSM we 
computed the hardness ratio map and we observed the same 
variation in the cluster centre, in the direction of the tem- 
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perature map elongation from south-east (SE) to north- 
west (NW). To assess the significance of the temperature 
variations found in the temperature map, we focused on the 
spectral fits in specific regions, see Fig. [6l The regions were 
chosen from temperature map. These results are shown in 
Table [3l The temperature obtained for the regions is in a 
good agreement with that from spectral wavelet analysis. 
More interesting results were obtained for the regions 9, 10 
and 3, where the temperatures are 7.81^0^93 keV, 12.0112^31 
keV and G.lGlJ'gy keV, respectively. The regions 9 and 3 are 
the temperature minima and the region 10 has the highest 
of temperature. The temperature peak does not coincide 
with the maximum of cluster emission. From our spectral 
analysis we established that the temperature variations are 
approximately 4 keV in the cluster center. It is important to 
understand whether these variations of temperature would 
give any significant effect on the total mass. 



Table 3. The spectral fits results in the chosen regions, 
with free abundances. Three EPIC cameras were used. 



region T{keV)/xted 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 



6.60_i.5o 

g^g7+0.87 



+3.14 /Q 77 
0.69/0-95 

6.16^0 87/0.69 

10.28+116/0.80 

9.63l?*^/1.04 

11.12tJ;^i/0.90 

13.77t4n/0-66 

9.25i;°Vl-01 
7.8llJ°i/0.77 
U.Oltitl/OM 




Fig. 6. The white contours illustrated the chosen regions 
from which spectra were extracted. The regions were de- 
fined on the basis of temperature map. The green contours 
are the chosen sectors, to extract the spectra for four tem- 
perature profiles. 



4-^ 



.10 1.00 

R, Mpc 



).01 0.10 1.00 

R, Mpc 



Fig. 7. The solid points are the total temperature profile. 
The dotted points are the obtained temperature profiles 
from sectors: a) the 'hot' temperature profiles ( the SE and 
NW directions) b) the 'cold' temperature profiles (the SW 
and NE directions) correspondingly. The 'cold' temperature 
profile (SW) in panel b) displays a local peak because of 
the offset of the temperature peak visible in Fig[6] 



We observed the increase in temperature profile com- 
pared with the total temperature profile in the directions 
of SE and NW and the decrease in SW and NE directions. 
The maximum temperature variations were observed in the 
cluster centre. We used these temperature profiles to deter- 
mine the variation on the total mass profiles. 



4.5. Temperature profiles in difFerent directions 

To determine the influence of these variations on the to- 
tal mass we extracted the temperature profiles in the four 
different directions. The SE and NW profiles show high 
temperatures whereas in the SW and NE direction they 
are lower. Fig. [6] shows the sectors and directions. The ref- 
erence point of sectors was chosen at the cluster emission 
centre. The spectra were extracted in each sector up to 4' 
(1.54 Mpc). The widths of the annuli were chosen similar 
to that of the total temperature profile. We obtained the 
temperature profiles in each direction. Fig. [7^ shows the to- 
tal temperature profile and two 'hot' temperature profiles. 
Fig. [7)3 shows the total and 'cold' temperature profiles. 



5. Mass analysis 

To calculate more precisely the total mass profile up to 
R200 we used the temperature and density profiles obtained 
above assuming hydrostatic equilibrium. To estimate the 
influence of the temperature variations on the total mass, 
we computed the total mass obtained with the different 
temperature proflles presented in Fig|5] and corresponding 
f3- parameters. So we performed a surface brightness anal- 
ysis in each chosen direction and we obtained /3-model for 
each proflles (see Table S]). 

Fig.[5]shows the obtained surface brightness proflles. We 
observed the different perturbation for each direction more 
signiflcantly in NW and SW direction from cluster centre 
to 500 kpc. 
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Fig. 8. The surface brightness profiles in different direc- 
tions. Solid line is the best fit of the total radial profile. 



Table 4. The /3-model fitting results in different directions 



direction 


/3-fit 


rc(Mpc) 




NE 


0.94 ± 0.01 


0.35 ± 0.01 


1.05 


SE 


0.65 ± 0.01 


0.18 ± 0.01 


1.22 


NW 


0.72 ± 0.01 


0.25 ± 0.01 


1.68 


sw 


0.77 ± 0.01 


0.31 ± 0.02 


2.3 



5.1. Cluster gas mass 

Using the electron density, we calculated the total gas 
density of the ICM and the ICM mass by integrating 
Eq.(3). We obtained a central electron density of rieo — 
(7.85 ± 0.01-)10~^cm~'^ assuming a temperature of 8.85 
keV and using the /3-model parameters obtained from the 
Nevalainen's background subtraction, given in Table[T] The 
total gas mass is shown in Table [5] 

5.2. Cluster total mass 

Assuming the spherical symmetry and the hydrostatic equi- 
librium we calculated the gravitational mass of the cluster 
CL0016+16. 

In the first step of the mass analysis we assumed that 
CL0016+16 is isothermal with the temperature of 8.85 keV 
and took into account error bars on the temperature profile, 
P and Rc parameters. The results are presented in Fig[51 

In the second step of the total mass determination we 
used obtained temperature profiles. Using the hydrostatic 
equilibrium approach, the /3 parameters obtained with the 
subtraction Nevalainen's background and the temperature 
profile we calculated the cluster total mass with: 

Mtot{< r) = ^jj—r' - m^^) (4) 

Gfirrip \ar r'^ + J 

The mass profile was calculated using the Monte Carlo 
method, which takes the obtained parameter for the gas 
density profile and the measured temperature profile as an 



Table 5. The results of CL0016+16 mass analysis using 
different temperature profiles and cluster parameters. 



parameter 


M(< i?200) 

(xlO'^'Afo) 


M(< Esoo) 
(x10"Mq) 


Mgaa 


2.24 ±0.06 


1.41±0.04 


Mtot,T=8.81 keV 
Mtoi,T(r) 


14.0±0.6 
11.9±2.0 


8.4±0.3 
8.0±1.0 


Mtot,T{r)NW hot 
Mtot,T{r)NE cold 

^ T(r\'^K tint 

Mtot,T{r)SW cold 


11.2±2.6 
11.1±2.4 
14.0±2.8 
12.0±2.5 


8.0±2.1 
8.1±1.9 

9.0±2.0 


VT 

Mtot,T = 8.85 keV 


14.42 




neo(xlO"^cm"^) 

f as 

LxIxlO-'^erg/s) 


7.85±0.01 
0.16±0.01 
5.1± 0.1 






1 1 

R, Mpc R, Mpc 



Fig. 9. The integrated total mass profile of CL0016-)-16. a) 
The bold solid line, showing the total mass profile, was ob- 
tained using the total temperature profile. The dot-dashed 
line is the total mass profile obtained using constant tem- 
perature, b) The bold solid line is the total mass profile 
obtained using the total temperature profile, c) The solid 
line is the total mass profile obtained using the 'cold' tem- 
perature profiles, d) The dotted line is the total mass profile 
obtained using the 'hot' temperature profiles. 



input. The method of [Neumann fc BohringeJ (jl997[ ) was 
used, which allows a transformation of the error bars of the 
temperature profile into error bars of the mass profile. 

Table [5] shows the results of our total mass determina- 
tion. Fig[5] shows the total mass profiles obtained using the 
constant temperature and assum ption temperature profiles . 

Using a similar approach of iKotov fc VikhlininI (|2005f ) 
we find the same result on the total mass. We extrapo - 
lated the temperature profile of iKotov fc VikhlininI (|2005() 
to i?200, obtaining total mass of (11.7± 1.7) • IO^^Mq, which 
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compare very well with the value of (11.9 ± 2.0) • 10^** Af© 
using our total temperature profile. 

In order to check the influence of temperature variations 
in the cluster centre on the total mass profiles, we deter- 
mined four temperature profiles in each direction and the 
total temperature profile. Using the equilibrium approach 
and the obtained /3 parameters the total mass was calcu- 
lated for each temperature profile. The values of the total 
mass for each temperature profile are shown in the Table [S] 
up to i?200 and i?5oo- These temperature variations does not 
significantly affect the total mass, AM < 20% at the i?2oo 
in comparison with the total temperature profile. The main 
contribution to the mass errors in these profiles came from 
the determination of temperature, which is less constrained 
in sectors than in the full annuli. It is noticeable anyway 
that the total mass estimates in sectors are all within 1 a 
errors. 

The total mass profiles for each temperature profile and 
total temperature profile are shown in Fig [HI We found that 
in the cluster's central parts the total mass profiles obtained 
with the 'cold' temperature profiles are lower than the to- 
tal mass profile obtained with the total temperature. At 
R200 we found the predominance of the total mass in the 
southern part of cluster. 

We calculated also the gas mass fraction of 0.16 ± 0.01 
which it is simply the ratio of the ICM mass to the total 
mass. 



5.3. Total mass from self similarity 

Also, the total mass was estimated from M ~ T scaling 
relation using the best fit of the mean temperature in 
5' radius at the cluster redshift, and in ACDM cosmol- 
ogy. Assuming structural similarity and the virial theorem, 
we provide a scaling relation between virial mass, radius 
and the overall X-ray temperature Tx: Mvt200 / Rvt200 oc 
Tx- This relation corresponds to a fixed density contrast 
at redshift z and is derived using the virial theorem: 
MyT2oo/(4/37rpc^^T2oo) = 200. This leads to the weh 
known scaling relation: 



VT200 



OC (1 + z)-3/2r3/2 



(5) 



To compute the scaling total mass Mvt2qo we used the 
n ormalization factor obtai ned from numerical simulations 
of iBrvan fc NormanI (jl998D . rescaled at Rvt200- The result 
for MvT200 is in the TableO From the virial theorem Eq.(6) 
and from the hydrostatic equilibrium Eq.(4) we obtained 
similar results on the total mass. 



6. Scaling properties 

6.1. Cluster profile 

To study the influence of the cluster properties on the self- 
similarity we compared the scaling and structural proper- 
ties of distant 'hot' galaxies cluster CL0016+16 with the 
scaling properties obtained from nearby galaxies clusters. 

The self similarity model is based on the simple as- 
sumptions for the cluster for mation, derived from th e top- 
hat spherical collapse model (jBrvan fc NormanI . |1998[ ) . The 
virialised part of a cluster present at a given redshift cor- 
responds to a fixed density contrast as compared to the 
critical density of the universe at that redshift. This model 
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Fig. 11. a) The dotted lines are the scaled emission mea- 
sure (ScEM) profiles of the nearby clusters. The radius is 
normalized to i?2oo- The bold solid line shows the ScEM 
profile of CL0016-f 16. b) The solid line is the mean disper- 
sion from ScEM profiles of the nearby clusters. The bold 
sohd line is the SEM profile of CL0016+16. 



makes the definitive predictions in terms of the evolution of 
cluster properties. We considered the scaled emission mea- 
sure profiles and the Lx — T relation. 

We used i?2oo for the normalization on the radius. The 
self-similarity of cluster implies that the scaled emission 
measure profile (ScEM) should be identical for all clusters: 



ScEM{r/R 



200) 



47r(l + z)45(r/i?20o) 



A(T,z)a 



:E{z)f 



(6) 



ScEM{r / R200) were calculated directly from the ob- 
served surface brightness distributions, S{r) is the obtained 
surface brightness profile. We applied Eq.(7) to all cluster 
profiles using the best fit temperature, corresponding red- 
shift, and the emissivity A{T,z). 

Using ACDM cosmology, we compared the ScEM pro- 
file of CL0016+16 observed with XMM-Newton in the 
0.3-2 keV band with the ScEM profiles of twelve nearby 
galaxy clusters observ ed with ROSAT in the 0.5-2 keV 
band (jNeumannl . l2005t ). The nearby galaxy clusters were 
obtained in the temperature range (1.7 < kT < 8.5) keV 
(but most of them in the 5-6 keV range). The relative error 
in the cahbra tion of XMM/EPIC and ROSAT/PSPC can 
be neglected (|Arnaud et al.l . [2003 ). FigfTT] shows the aver- 
age ScEM profiles of nearby galaxy clusters and the ScEM 
profile from CL0016-I-16. 

We confirm that, in the cluster centre, the ScEM profiles 
show a large dispersion, which are commenly be explained 
by the non-gravitational processes. The ScEM profile of 
CL0016-f 16 is lower in the centre in comparison with aver- 
age ScEM profile, which is the consequence of the absence 
of a cooling core. In outer regions (r > 0.2i?2oo) the ScEM 
profile of CL0016+16 is higher than the average ScEM pro- 
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Fig. 10. a) The X-ray emission from the XMM-Newton data b) The X-ray residual from the XMM-Newton data in the 
cluster centre, was obtained with the help of 2D /3-model and c) The obtained X-ray residual from the Chandra data. 



file of nearby galaxy clusters. We compared the total emis- 
sion from CL0016-t-16 and from nearby galaxy cluster at 
the radii 0.2i?2oo > r > i?200- The total scaled emission 
of CL0016-I-16 is by a factor of 1.6 stronger than the total 
emission from nearby galaxy clusters. 

Finally, these results can be explained by cluster prop- 
erties of CL0016-I-16, in particular by the absence of cooling 
core and high luminosity which in merge activity. 

6.2. Lx-T relation 

We tested the scaling Lx — T relation for distant lumi- 
nous galaxy cluster with possibility of merger. We com- 
puted the bolometric luminosity Lx within the i?2oo (see 
Table [5]). The observed count rate in the 0.3-4.5 keV band 
was obtained. The integrated surface brightness profile was 
converted to bolometric luminosity using the best fit 
model from Nevalainen's background subtraction and the 
instrumental response. The error in Lx includes the statis- 
tical errors on the count rate and te mperature. We tested 
the Lx — T using scaling relation from lArnaud and EvrardI 
(fT999l ). Lx oc This relation was obtained for sam- 

ple clusters with weak or absent cooling flow signatures, 
so it can be used for CL0016-I- 16. Taking the simple mean 
temperature, scaling relation bv lArnaud and EvrardI ()1999[ ) 
for our cosmology, and normalization evolution factor to 
our distant cluster. We found that CL0016-I-16 is more lu- 
minous by a factor of 2.1 for this temperature. These re- 
sults confirmed the possibility of merger in maximum core 
collapse for CL0016-I -16, because the numerical simulation 
([Randall et al.l . l2002f) suggest that the effect of the merger 
is to boost the luminosity and the temperature for a short 
time and this can have an obvious effect on scaling relation. 

6.3. 2D (3-analysis, X-ray residuals 

To better understand the status of the cluster dynamics 
in the centre we performed a two-dimensional fit to the 
cluster surface brightness, using a modified /?-model that 
allows for two different core radii along the two principal 
axes of cluster image ellipse. We fited the surface brightness 
distribution of CL0016+16 obtained from XMM-Newton 



data with a 2D /3-model and we quantified the deviation 
from this model. 

For the XMM-Newton data we used three cameras, with 
the modelling and subtraction of background from image 
(same method were used in section 3.1.3) and with point 
sources and gap corrections. 

To confirm XMM/Epic 2D (3 analysis, we used also the 
Chandra data to compare the residuals at better spatial 
resolution. We are confident of the reality of these resid- 
uals since we found very similar results with XMM/EPIC 
and Chandra/ACIS. Reduced cleaned ACIS event list was 
downloaded directly from Chandra X-ray Center (CXC). 
The analysis was performed with CIAO software. The im- 
ages were smoothed with Gaussian filter before fitting. We 
fitted 2D /3-model in the region 2' (0.8 Mpc « i?5oo). The 
best-fit parameters are listed in Table [51 AH parameters 
were assumed free in the fit. 

The residuals are shown in the Fig.lTOl From the XMM- 
Newton and Chandra data we obtained very similar results. 
But we obtained a different core radius because we used a 
different satellite, with a different PSF. Notice, that we do 
not take the PSF into account, this explain the different in 
the core radii. Two very similar residual maps established 
the perturbations at the centre over at least 1' radius ( 0.4 
Mpc). In both case the perturbations in the cluster cen- 
tre were observed with 20% deviation from the maximum 
cluster emission. 

After radial projection, we obtained a rather good (3 fit 
for total surface brightness profile, but we observe pertur- 
bations in 2D. At first glance, CL0016-I-16 looks rather ho- 
mogeneous in brightness and not in temperature. We now, 
after the 2D f3 subtraction, enhance the brightness struc- 
tures in the inner core. We do observed perturbations in 
the temperature and in the surface brightness distributions 
in the cluster centre. Quantitatively, these 2D f3- model 
residuals do not lead to strong variations in density profiles 
(thus in the mass estimate in the HE scheme) , but they are 
of importance to qualify the relaxation status of the inner 
core of CL0016-I-16 and thus the limit of the hydrostatic 
equilibrium hypothesis. 

If CL0016 is at the maximum core collapse phase, it 
will present very similar residuals as observed in numerical 
simulations around this phase. See, for example, figure 4 
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in iRicker and Sarazinl (|200lD , where they present a 1 to 1 
mass ratio collision with high impact parameter. Following 
that scheme, we can estimate the time since core collapse 
passage by more precise comparison with these simulations. 
CL0016-I-16 is 2 time too luminous as what the simulation 
found for a cluster at 0. 01 Gyr around its maxi mum core 
collapse (see figure 5 in (jRicker and Sarazinl . l200in ). 

Table 6. The best-fitting results of the 2D /3-model, from 
the XMM-Newton and Chandra data 



Parameter 


Best-fitting 
XMM 


Chandra 


i?ci(Mpc) 


0.291 


0.237 


Rc2{Mpc) 


0.364 


0.295 


P 


0.79 


0.79 


Pa 


2.24 


2.33 


RA 


00:18:33 


00:18:33 


Dec 


16:26:06 


16:26:11 



7. Discussion 

In previous works, CL0016-I-16 was always considered as 
a very massive relaxed cluster. In our present analysis, we 
found the strong evidence that, at least in the centre, the 
cluster is not relaxed. From the weakest to the strongest 
argument favouring the merger scenario, we found: 

— CL0016-f 16 was about a factor of 2 too bright, when 
compared to the expected luminosity of the Lx — T 
relation. 

— We detected similar significant residuals after subtrac- 
tion of a 2D /3-model fit using both XMM-Newton and 
Chandra observations, especially in the inner region 
(r < 1' « 400 kpc). Notice that the X-ray maximum 
was not onto the centre of the D /3-model fit. 

— The temperature map of the 2.5' region (i.e. up to 1 
Mpc) clearly looks like an equal mass merger at first 
maximum core collapse showing significant azimuthal 
variations. 

In addition, CL0016-I-16 is also well known to host a 
strong r adio halo which also argues in favour of a merging 
cluster ([Giovannini fc Ferettil . |2000D . Also we obtained a 
high velocity distribution of galaxies for this cluster. Using 
NED data, we obtained a <7„ = 1800 km s~^ for the 150 
galaxies at CL0016+16 redshift. The maximum of this map 
is offset by 0.4' (141 kpc) with the X-ray one. We tried to 
quantify the impact this could have onto the total mass es- 
timation and found that depending on the adopted temper- 
ature profile, (but always in the framework of Hydrostatic 
Equilibrium), the total mass may vary by something like 
20%, which is much bigger than the quoted error which 
are around 5% (the errors in the temperature profile is the 
main contributor to this error budget). 

The main results of our detailed analysis of the galaxy 
cluster CL0016-I-16 from the XMM-Newton data can be 
summarized as follows. 



We performed a detailed image and spectral analysis of 
the XMM-Newton data using three different backgrounds. 
In all cases we detected the cluster emission up to i?2oo 
with the XMM-Newton data. Also using the template of 
simulated cluster we checked the possibility of detection 
with XMM-Newton data up to i?2oo (see Appendix E]). All 
obtained surface brightness profiles give a good fit with the 
/3-model. 

We studied the dynamics of CL0016+16 from the de- 
tailed spectral analysis. The global temperature estimated 
in 5' is 8.81 keV, but the spectral study shows that this 
cluster is not isothermal. We obtained the temperature pro- 
file which decreases in the outer regi ons, our results are 
i n a g ood agreement with results from lKotov fc VikhlinirJ 
(|2005f ). The temperature map shows the asymmetry in the 
radius 1', the temperature maximum (T = 12.0 ± 1.5 keV) 
is not superposed on the centre X-ray cluster emission and 
is located to south-western direction; also we observed the 
cold regions (T = 6.2± 1.2 keV) to north-east. To study the 
dynamics of CL0016-I-16 more precisely, we chose several 
regions using the temperature map. We extracted temper- 
ature in the regions, found no spherical symmetry in tem- 
perature and the presence of the 'cold' and 'hot' regions in 
the cluster centre. 

We tested the influence of the temperature variations on 
the cluster total mass using the obtained 'cold' and 'hot' 
temperature profiles. We calculated total mass within i?200 
and i?5oo, the temperature variations do not significantly 
impact on the total mass, AM < 20% at i?200- But the 
temperature variations in the centre influence mainly on the 
inner part of total mass profile. To better understand the 
possible dynamics in the cluster centre we calculated 2D /3- 
model from XMM-Newton and Chandra observations. The 
same residuals were obtained in the cluster centre, which 
may be an indication of merging. 

We tested CL0016-I-16 on the self similarity theory. 
ScEM of distant galaxy cluster CL0016-I-16 is higher than 
the mean ScEM profile of nearby cluster in external re- 
gions. In the centre we saw larger dispersion, which can be 
explained by non-gravitational process es. We tested Lx—T 
relation using scaling relation from lArnaud and EvrardI 
()1999D and concluded that CL0016-I-16 is more luminous 
by factor of 2 for this temperature. It is also the argument 
for maximum core collapse in CL0016+16. 
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changing the radial profile (below, we refer to this cluster as 
simulated cluster). To obtain the cluster emission from the 
simulation, we created the photons using the Monte Carlo 
approach. We calculated emission measure per cell using 
the cell size, density and the solid angle w hich depends on 
the used cosmology (jBourdin et al .1 . 1200^ . 

The photons were obtained from cosmology simulation 
for the simulated cluster. We convolved these photons with 
XMM-Newton response. In our analysis we took into ac- 
count the PSF effect, vignetting correction function, mask 
of CCD cameras (gap) and added the randomly selected 
events from background data of A. Read. The same data 
reduction was performed as for CL0016-I-16. Fig. 4 shows the 
resulting surface brightness profile for the simulated cluster. 
Notice that Cluster 6 is more cuspy than CL0016-I-16, its 
/3-model parameters are /3 = 0.99 and the core radius rc — 
0.24 Mpc. We still detect the emission of simulated cluster 
up to i?2oo but with less significance, see Fig lA.ll This is a 
rather robust confirmation that XMM-Newton can detect 
the bright cluster emission up to i?200- 




Appendix A: Possibility of detection up to i?2oo 
with XMIN/I-Newton 

In order to check the possibility of the cluster emission de- 
tection up to i?200 with XMM-Newton data we decided to 
treat a template of similar cluster fro m the Hydro N- body 
cosmology simulation, code RAMSES (jTevssierj . [2003 ). 

The leading idea to use numerical simulation was only 
to check the ability of XMM/EPIC to detect such a massive 
cluster at this rcdshift up to its virial radius taking as pre- 
cisely as possible all the observational effects into account 
in particular, the relationship between the background, the 
instrumental response and the hot temperature of this clus- 
ter. The scaling applied onto Cluster6 directly follow the 
M-T relation and allow to properly study the instrumental 
effect and limits. 

We want to obtain similar cluster using the adaptive 
mesh refinement 3-D hydrod ynamical cos mology simula- 
tion of structure formation ( Tevssierl . l2002f ). Cluster 6 was 
formed at a temperature of 3.4 keV. 

Using the self-similarity theory from Cluster 6 we scaled 
its temperature and mass to those of CL0016-I-16 without 



Fig. A.l. The surface brightness profile for similar cluster 
of CL0016-I-16 obtained from the Hydro N-body cosmology 
simulations. 



